小半径弯管弯胀复合成形研究

为解决小半径弯管内侧起皱、外侧减薄的成形缺陷,提出了一种弯胀复合成形的方法。通过几何关系解析了初始弯曲半径和初始管坯直径与最终管径和最终弯曲半径之间的关系,确定了初始弯曲半径的取值范围。采用有限元与实验相结合的方式研究了不同初始弯曲半径对零件壁厚分布的影响,分析了胀形过程中零件弯曲处及直壁处的变形规律。最终发现:初始弯曲半径越大,零件的最终减薄越小,胀形过程中弯曲处内外侧发生均匀减薄,直壁处内侧减薄明显大于外侧减薄。     

多种工质下螺旋波等离子体源波场结构数值模拟

针对螺旋波等离子体放电机理，开展了多种工质条件下的螺旋波放电等离子体内波场结构数值模拟研究。计算发现:氦气等离子体的Er分量在径向边界处的峰值更为突出，有利于等离子体在径向的输运，波电场径向分量决定了电流密度径向分量在内部的表现。在0.266 Pa和1.064 Pa两种气体压强条件下，通过波场结构验证了气压对于波阻尼影响的结论。波场结构是螺旋波在等离子体内传播以及能量沉积的微观体现，研究螺旋波波场结构是揭示其高电离效率的重要途径。初步探索了功率耦合机制，为实验系统优化及实验方案设定奠定理论基础。     

Fast spacecraft solar radiation pressure modeling by ray tracing on graphics processing unit

A novel method is presented to evaluate on the graphics processing unit (GPU) the force and torque on a spacecraft due to solar radiation pressure. The method employs efficient ray tracing techniques, developed in the graphics rendering discipline, to resolve spacecraft self-shadowing and reflections at faster than real-time computation speed. The primary algorithmic components of the ray tracing process which contribute to the method’s computational efficiency are described. These components include two-level bounding volume hierarchy acceleration data structures, fast ray to bounding box intersection testing using the slab intersection algorithm and fast triangle intersection testing using the Möller-Trumbore algorithm. Spacecraft material optical properties are represented as a combination of Lambertian diffuse and ideal specular reflections. Both diffuse and specular ray-surface interactions are modeled. The approach is implemented using C++ and OpenCL and executed on a consumer grade GPU. Model validation is presented comparing ray traced force and torque values to the same quantities produce by a faceted analytic model. Numerical results illustrate the impact of self-shadowing on the force and torque calculation, and demonstrate the fast computational speed that is enabled with this implementation.     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

回流燃烧室流动特性试验

为了揭示有/无燃烧状态下燃烧室热态和冷态流场的特征和流动特性，针对某型回流燃烧室单头部试验件，使用粒子图像测速仪（PIV），测量燃烧室燃烧状态下不同截面处的热态流场，以及没有燃烧状态下不同截面处的冷态流场，探讨不同总压损失系数对回流燃烧室热态/冷态流场特征及流动特性的影响。研究表明:随着总压损失系数的增大，冷态条件下各截面流场结构基本保持不变，如射流孔穿透深度、射流角度、回流区位置及大小、流线等基本保持一致，但是各位置点速度大小逐渐增大。热态条件下各截面流场随着总压损失系数增大，流场结构也基本保持不变；相同总压损失系数时，热态流场与冷态流场存在差异，燃油喷射与气流的相对运动将会对燃烧室头部的流场结构造成影响，速度较冷态流动时略微增大。      

弹片热机械疲劳试验方法研究及设计

针对弹片的热机械疲劳（TMF）试验要求，采用机械设计技术、机电技术、冷却技术、计算机技术和数据采集技术，提出了压力载荷和热载荷的加载方法，建立了2种载荷5个试验件的并联试验控制系统，并设计了弹片热机械疲劳试验器。试验结果表明:该系统能够同时模拟服役环境下弹片的压力载荷和热载荷。利用该试验系统进行了弹片的热机械疲劳试验，试验结果再现了弹片在服役状态下的失效模式。试验系统具有良好的重复度、较高的加载频率和加载精度。压力载荷最大相对误差为2.4%，绝对误差小于0.5N。温度载荷最大相对误差为3.55%，最大绝对误差为3.89℃。      

基于两相流模型的挤压油膜阻尼器空化流场特性数值模拟

基于气液两相流理论中的Mixture模型和Zwart-Gerber-Belamri空化模型,采用动网格技术建立了两端开口中心槽供油型挤压油膜阻尼器三维非定常空化流场求解模型。数值模拟表明:随着阻尼器内环进动,考虑两相流动的油膜低压区在中心槽的两侧产生两道对称的条状负压带,条状带内具有较高的气相体积分数;阻尼器低压区的压力和气相体积分数对进油孔位置十分敏感,阻尼器油膜力和流场气穴比的变化频率与进油孔数密切相关。内环同心进动半径、进动频率对阻尼器空化流场影响的数值计算表明,进动半径和频率的增大均会使得流场内的空化现象加剧,同时气穴比的相位滞后现象愈加显著。      

风-雨耦合环境超高层三塔连体建筑雨压分布研究

以中国东南沿海某超高层三塔连体建筑为研究对象，以风-雨双向耦合算法为核心，基于计算流体动力学技术采用连续相和离散相模型进行风场和雨场的迭代模拟。首先基于9种风雨组合工况进行三塔连体建筑非定常脉动风场模拟，探讨超高层三塔连体建筑平均风压分布、表面速度流线和流场干扰机理。然后对比研究不同风雨组合工况下主塔表面雨滴附着数量、雨滴冲击力和雨压系数的分布规律，揭示风-雨耦合场中结构表面速度流线、雨滴运行轨迹和最终速度的作用机理。最后提炼出超高层三塔连体建筑最不利风-雨组合工况，并给出对应的雨压系数取值建议。研究表明：风-雨耦合环境下超高层连体建筑迎风面雨荷载作用最为显著，此时雨荷载与风荷载最大比值可达23.81%，局部测点最大雨压系数达到0.301，100 a重现期风速和强大暴雨组合为风-雨耦合作用的最不利组合工况。     

带轴对称台阶的圆锥高超声速边界层转捩试验

在Ma 6风洞内，通过高频脉动压力测试技术和基于纳米粒子示踪的平面激光散射(NPLS)技术，分别对带前向、后向轴对称台阶的圆锥高超声速边界层转捩进行了试验研究。采用功率谱密度分析和互相关计算等方法对脉动压力数据进行分析，得到了边界层中扰动波的发展规律，定量分析了第二模态波的相关参数。结果显示：两种模型中第二模态波在沿流向向下游发展的过程中，其幅值均先增大再衰减、特征频率均逐渐减小；特征频率和传播速度整体上均随雷诺数的增大而增加(后台阶模型中特征频率由100 kHz增至196 kHz，前台阶中则由 97 kHz 增至174 kHz)、波长变化规律则与之相反(后台阶中由6.35 mm降至4.54 mm，前台阶由7.35 mm降至 4.66 mm )；后台阶模型中第二模态波初次出现位置比前台阶中更靠近上游，边界层转捩位置较前台阶前移。将NPLS结果与高频脉动压力测试结果进行对比，两者吻合较好。     

Numerical investigation on flow and heat transfer characteristics of supercritical RP-3 in inclined pipe

Numerical simulations of flow and heat transfer to supercritical RP-3 through the inclined tubes have been performed using LS k–e model embedded in Fluent. The physical properties of RP-3 were obtained using the generalized corresponding state laws based on the fourcomponent surrogate model. Mass flow rate is 0.3 g/s, system pressure is 3 MPa, inlet temperature is 373 K. Inclination of the inclined pipe varied from -90° to 90°, with heat flux varied from 300 k W/m~2 to 400 kW/m~2. Comparison between the calculated result and the experimental data indicates the range of error reasonable. The results of ±45° show that temperature inhomogeneity in inclined pipe produce the secondary flow in its cross section due to the buoyancy force. Depending on the strength of the temperature inhomogeneity, there will be two different forms of secondary flow and both contribute to the convective heat transfer in the pipe. The secondary flow intensity decreases when the inhomogeneity alleviates and thermal acceleration will play a leading role. It will have a greater impact on the turbulent flow to affect the convective heat transfer in the pipe. When changing the inclination, it affects the magnitude of the buoyant component in flow direction. The angle increases, the buoyancy component decreases. And the peak temperature of wall dominated by the secondary flow will move forward and increase in height.